arXiV: astro-ph/0112027v 1 3 Dec 2001 


The Astrophysical Journal, ???:???-???, astro-ph/ 0 112027 
Preprint typeset using L?T^ style emulateapj v. 25/04/01 


THE DENSITY AND LOCATION OE THE X-RAY ABSORBING GAS IN NGC 3516 

Hagai Netzer,^ Doron Chelouche,' I.M. George,TJ. Turner,D.M. Crenshaw,"^ S.B. 

Kraemer,^ K. Nandra,^ 

Accepted by The Astrophysical Journal 

ABSTRACT 

A new Chandra observation and archival observations by ASCA are used to investigate spectral variations in 
the Seyfert 1 galaxy NGC 3516 over a period of 7 years. A large change in flux (factor ~ 50 at 1 keV) is observed 
between an ASCA observation in 1994 and the Chandra observation in 2000, with the source close to the all-time 
maximum and minimum X-ray flux states, respectively. We find the variations in the observed flux and spectra at 
these epochs to be consistent with a constant column density of line-of-sight material reacting to changes in the 
ionizing continuum. The data from the two epochs are consistent with a simple decrease (by a factor 8-10) in the 
luminosity of a constant 0.5-50 keV slope source and a line-of-sight absorber with an equivalent hydrogen column 
density of 10^^ ® cm^^. Intermediate luminosities, sampled during other A5CA observations, are all fitted by the 
same model with a very small change in spectral index (well below AT = 0.2). In addition, analysis of the long 
(360 ks) ASCA observation in 1998 shows clear “color” variations that are entirely consistent with this model and 
are interpreted as due to changes in the opacity of the absorbing gas. The data allow us to put a conservative upper 
limit of 60 ks on the recombination time which translates to a lower limit of about 2.4 x 10 ® cm ^ on the density 
of the recombining gas and an upper limit of about 6 x lo^^h-i cm on its distance from the central source. These 
are the best limits obtained so far on the density and location of the X-ray absorbing gas in a type-1 Active Galactic 
Nucleus (AGN). They indicate that the absorbing gas is different, in terms of its density and location, from the 
ionized gas commonly observed in type-II AGN. The Chandra ACIS/LETGS data also reveals a strong (EW=290 
eV), unresolved 6.4 keV iron line, a strong O VII 0.561 keV line and a marginally detected N VI 0.419 keV line. 
The former is interpreted as originating in a large column of gas of lower state of ionization seen in “reflection”, 
and is consistent with the spectrum at high energies at all epochs. The two others emission lines are probably 
emitted by the gas also responsible for the line-of-sight absorption. 

Subject headings: galaxies: active — galaxies: individual (NGC 3516) — galaxies: nuclei — galaxies: Seyfert 
— techniques: spectroscopic — X-rays: galaxies 


1. INTRODUCTION 

Most low-luminosity Active Galactic Nuclei (AGN) 
contain significant amounts of circumnuclear gas which mod¬ 
ifies the intrinsic X-ray spectrum by introducing strong ab¬ 
sorption features, around 1 keV. This gas, which has been 
named “Warm Absorber” or the “highly ionized gas” (hereafter 
HIG), must be photoionized and of moderately low temperature 
(- 10®K), as deduced from many detailed studies. Works by 
Reynolds, (1997), George et al. (1998, hereafter G98) and oth¬ 
ers show it to be present in some 50-70% of all low luminosity 
type-I AGN (i.e. broad emission line Seyferts). HIG systems 
appear to be less common in higher luminosity type-I objects 
but the statistics are less secure (George et al. 2000). Previous 
attempts to analyze the gas properties focused on the absorbing 
columns and the level of ionization. The last property is well 
described by the X-ray (0.1-10 keV) ionization parameter, Ux 
(Netzer 1996) which is the ratio of the 0.1-10 keV photon den¬ 
sity to the gas density. G98 showed that while objects contain¬ 
ing X-ray absorbing systems vary a lot in the column density 
of ionized gas, there is a surprising uniformity in Ux with most 
objects clustering around Ux=QA. 

As yet, there is no real indication of the location of X- 
ray absorbers. Plausible arguments have been made both for 


kpc scale and for much smaller distances (e.g. Krolik & Kriss 
1995, 2001; Netzer 2001 and references therein). Because of 
that, there is also no information about the mass of the emit¬ 
ting material and its kinetic energy. The best constraints can be 
obtained by high resolution X-ray spectroscopy yet, the fun¬ 
damental issues are not resolved in the handful of available 
sources (e.g. NGC 5548, see Kaastra et al. 2000; NGC 3783, 
see Kaspi et al. 2001; MCG-6-30-15, see Branduardi et al. 
2000 and Lee et al. 2001) despite of the wealth of emission 
and absorption lines discovered by the new Chandra and XMM- 
Newton observations. 

The target of our study, NGC 3516 (z=0.009) is at a dis¬ 
tance of 35.1hj^ Mpc where /175 = Ho/15 km s^^ Mpc^* (Per¬ 
mit et al. 1998). This is a moderate luminosity (L 2 _iokeV — 
jq42.5-43.0;j-^2 ^-1 1994_1998) type-I AGN with 

typical X-ray properties and a small column (Nf/f ~ 3 x 10^*^ 
cm^^) of Galactic absorption. The ASCA based studies of 
Kriss et al. (1996b), Reynolds (1997) and G98 indicated the 
HIG has a large absorbing column (Nh ^ 10^^ cm^^) and av¬ 
erage ionization parameter. Detailed observations with RXTE 
reveal large variability on long ('^ 100 days) and short (days) 
time-scales (Edelson et al. 1999, Markowitz and Edelson 
2001). Nandra et al. (1999) reported the presence of a rela- 
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tivistically broadened, variable, Ee 6.4 keV line. A recent pa¬ 
per by Guainazzi et al. (2001), based on BeppoSAX observa¬ 
tions, attributed the X-ray variations during 1996-1997 to large 
columns of cold material moving in and out of the line of sight. 

NGC 3516 has been extensively observed in most other 
wavelength bands. Most relevant to our study are the UV obser¬ 
vations (Walters et al. 1990; Kriss et al. 1996a, and references 
therein; Crenshaw et al. 1999; Huchings et al. 2001) showing 
strong, multi-component absorption lines superimposed on the 
broad emission lines like C IV A.1549 and N V A.1240. The ab¬ 
sorption lines are highly variable on months to years time scales 
and suggest gas outflow from the nucleus at a large range of 
velocities, 0-1300 kms^^ Walters et al. (1990), Kriss et al. 
(1996a) and Mathur, Wilkes and Aldcroft (1997) discussed var¬ 
ious aspects of this UV-X-ray relationship. 

This paper reports on new. Cycle 1 Chandra observa¬ 
tions of NGC 3516. The source happened to be close to its all- 
time X-ray minimum and hence the signal-to-noise (S/N) is not 
adequate for a detailed study of emission and absorption lines. 
However, the very low luminosity provided a unique opportu¬ 
nity to study the time-dependent properties of the absorbing gas 
and to obtain the first firm limits on its density and location. 
The new observations are described in §2, the physical model 
and the time-dependent properties are presented and discussed 
in §3 and a general discussion is given in §4. 


2. OBSERVATIONS AND REDUCTION 

Chandra observed NGC 3516 on 2000 October 6 for 
47 ks with the Low Energy Transmission Grating Spectrome¬ 
ter (LETGS) in the optical path. The “S-array” of the Advanced 
CCD Imaging Spectrometer (ACIS) was used as the focal-plane 
detector. A sub-array read-out was employed with a 1.1 sec 
frame time. Standard (ClAO v2.1) pipeline spectral extraction 
was performed on the data using the corresponding calibration 
files (CALDB v2.6). The data were screened for known bad pix¬ 
els and columns and “bad grade” events were discarded. The 
screening resulted in a net exposure of 44 ks. In order to cor¬ 
roborate our results from the LETGS ± 1 -orders, and to improve 
our signal-to-noise ratio at energies below 1 keV, we have also 
extracted and analyzed the zeroth-order spectrum. The total 
number of counts is 6800 for the zeroth-order and 8300 for the 
combined ±1 orders. Background is negligible for the grat¬ 
ing spectrum and was not subtracted to avoid corrupting the 
low S/N low-energy spectrum. Background has been subtracted 
from the zeroth-order spectrum. 

Pile-up of the zeroth-order spectrum, at low energies, 
may become important, despite the extremely low flux level of 
the source at the time of observations. Timing analysis, based 
on a 500 s binning of the zeroth-order data, shows a constant 
count-rate of 0.16 counts per 1.1 sec frame throughout the 47 
ks observation . Thus uncertainties due to pileup are well below 
the other observational uncertainties. To further test the flux 
calibration accuracy, we have made a very detailed comparison 
of the zeroth-order spectrum and the heavily binned E'-order 
spectrum. The two are shown in Eig. 1. The comparison shows 
a good agreement and complete consistency within the observa¬ 
tional uncertainties except for a small discrepancy over the 1-2 
keV range. This discrepancy has already been noted by others 
(e.g. Turner et al. 2001) and, despite considerable improvement 
in instrument calibration, is most probably the result of residual 
calibration problems with the E'-order data. 


3. THE CENTRAL CONTINUUM AND THE ABSORBING GAS 

NGC 3516 has been observed in several optical and IR 
wavebands (Pogge 1989; Miyaji et al. 1992; Permit, Wilson 
and Mulchaey 1998 and references therein). The central region 
show extended NLR emission, on a ~ 2 — 10 arcsec scale, with 
clear spiral-arms type structure. X-ray imaging by RG&4r-HRI 
has been reported in Morse et al. (1995) showing a marginally 
significant extended emission in the same general direction as 
the optical line emission. Our zeroth-order Chandra observa¬ 
tions clearly show signs of extended nuclear emission. This is 
beyond the scope of the present paper, which focuses on the 
new spectroscopic results, and is addressed in a separate publi¬ 
cation (George et al. 2002). 


3.1. Spectral analysis of the 1994 and 2000 observations 

Inspection of the new flux-calibrated data (Pig. 2) shows 
an unusual X-ray spectrum, with extremely hard continuum 
for a type-I source and a highly curved low energy spectrum. 
There are strong and broad absorption features at 0.6-0.9 keV, 
two strong unresolved emission lines (O VII 0.561 keV and 
Pe 6.4 keV) and a weak, noisy emission feature at 0.419 keV 
identified as N VI 0.419 keV. There is also a hint of a weak emis¬ 
sion line near 905 eV that could be due to Ne IX. While the 
latter raises the interesting possibility of a highly-ionized com¬ 
ponent, the S/N of this feature is so poor that we neglect it 
in our analysis. Table 1 gives a summary of the more impor¬ 
tant results of the new observations with their lo error bars. 



Fig. I.— A comparison of the zeroth-order (red points) and the P*-order 
spectra of NGC 3516. The combined ±I orders data were binned to a constant 
EjIsE = 30 resolution. The difference between instruments over the 1-2 keV 
range is likely due to the known LETGS calibration problem. The insert shows 
the 0.45-0.6 keV range with the large EW O VII line at 561 eV. 


Table 1: New Chandra results 

Flux(2-10 keV) 

1.35x10 “ ergs/cm'^/sec 

Flux(0.5-10 keV) 

1.55 X 10“*' ergs/cm^/sec 

EW(0 VII 0.561 keV) 

14±5 eV 

EW(Fe 6.4 keV) 

290±50eV 

Ff (continuum at 6.4 keV) 

(1.6 ±0.2) X 10“"^ photon/ cm^/sec/keV 

Flux(2-10 keV, 1994) 

7.35x10“" ergs/cm"/sec 
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Fig. 2.— ACIS/LETGS spectrum and model fit. The two first-order spectra 
are combined to an even resolution of E /AE = 50 and the model is rebinned 
to a similar resolution. The reflection component, with the assumed normal¬ 
ization of 32% of the total continuum flux at 6.4 keV, is shown in blue. The 
deviation between model and data over the 1-2 keV range is attributed to the 
known LETGS calibration problem (see Fig. 1). 

A comparison of the new observations with the 1994- 
April ASCA spectrum suggests a flux drop by a factor 6.4±0.5 
at 5 keV and a factor ^ 50 at 0.7 keV. Thus the spectrum has 
become harder than almost any other known type-I AGN spec¬ 
tra. These changes suggest either extremely large intrinsic slope 
variations or a large wavelength dependent change in the opac¬ 
ity of the line-of-sight gas. As shown below, the first of these 
possibilities seems to be inconsistent with the available obser¬ 
vations. The following is a detailed discussion of the second 
possibility based on detailed modeling of the central continuum 
and the absorbing gas. 

We have performed detailed photoionization calcula¬ 
tions in an attempt to explain the spectra of NGC 3516 at all 
epochs. Our model calculations use ION2001, the 2001 ver¬ 
sion of the photoionization code ION (Netzer 1996; Kaspi et 
al. 2001 and references therein). The code includes a self- 
consistent treatment of the ionization and thermal structure of a 
steady-state ionized gas exposed to an external radiation source. 
Atomic data include the most recent cross sections for the more 
abundant elements and all of their strong lines. They also in¬ 
clude new f-values for a large number of iron L-shell lines (see 
Kaspi et al. 2001), newly calculated ls-2p transitions (Behar 
and Netzer 2002) and dozens of M-shell iron transitions (Behar, 
Sako & Kahn 2001). These models are more sophisticated than 
any other photoionization models applied to this source because 
of the large number of lines, the inclusion of both emission and 
absorption features and the accuracy of the atomic data. 

Throughout this work we used models of constant den¬ 
sity («//) line-of-sight clouds that are specified by the contin¬ 
uum spectral energy distribution (SED), the equivalent hyd- 
grogen column density of the gas (Nh) and the “oxygen ion¬ 
ization parameter” {Uoxi defined over the 0.538-10 keV range, 
see Netzer 2001). As explained below, fitting of the 4-8 keV 
continuum of all X-ray observations (an energy range which 
is not affected by absorption) requires two somewhat differ¬ 
ent SEDs with photon indices of r(0.5 —50 keV) = 1.7 and 
r(0.5 —50 keV) = 1.9. The assumed SEDs are specified by 
broken power-laws between the energies of 0.01, 0.1, 0.5 and 
50 eV. In both cases r(0.01-0.1 keV)=2. Eor the r(0.5 - 


50 keV) = 1.7 SED case, we have r(0.1-0.5 keV)=3.3 and for 
the r(0.5 -50keV) = 1.9 SED, we need r(0.1-0.5 keV)=3.0. 
The poor S/N below 0.5 keV does not allow clear determina¬ 
tion of slope (or SED shape) at those energies and the spectral 
indices over those ranges represent, for example, the smooth 
exponential distribution of an EUV source. Eor clarity and ref¬ 
erences to other commonly used definitions of the ionization 
parameter we note that for the r(0.5 —50 keV) = 1.9 case, 
17(13.6 eV):t/^;t/o,,=178:13.1:l and for the r(0.5-50keV) = 
1.7 SED, t/(13.6 eV);t/^;t/oj=240:15.4:l, where Ux is defined 
over the 0.1-10 keV range and t/(13.6 eV) over the hydrogen 
Lyman continuum. 

The overall normalizations of the above SEDs are based 
on two observed points, at 1470A (Kraemer et al., 2002) taken 
on 2000 October 1, and near 7 keV (our own observations) 
taken 5 days later. We assume that this UV-X-ray flux ratio 
has not changed throughout the 7-year period and defer the 
more detailed discussion to Kraemer et al. (2002). In the X-ray 
band, the model results are not very sensitive to this constant- 
SED assumption since the ionization of the X-ray absorbing 
gas is insensitive to the UV radiation. While the chosen 7 
keV point is not affected by absorption, this is not the case 
for the 1470A point which is affected by reddening. This is¬ 
sue has been discussed by Kriss et al. (1996a) in their analysis 
of the 1995 observations of the source. Our adopted value cor¬ 
responds to Eb-v = 0.1 mag., and is somewhat larger than the 
value adopted by Kriss et al.. It is based on the unusually curved 
optical-UV continuum of NGC 3516 as observed by HST in 
1998 (see Edelson et al. 2000) and may, in fact, represent a 
lower limit to the line-of-sight reddening. A more detailed dis¬ 
cussion of the SED and the issue of reddening is deferred to a 
future publication (Kraemer et al. 2002). 

The following analysis is based on fitting the ASCA and 
Chandra spectra using standard minimization. Throughout 
this part we adopt the acceptability criteria suggested by G98 
i.e. a model is considered to be an adequate description of the 
data if P(x^ |dof) 0.95. We note that two of the model param¬ 
eters, the covering fraction of the emitting gas and the EWHM 
of the absorption lines, cannot be constrained by the data yet 
they can affect the calculated 'z}. After testing several cases, 
we have fixed those to be a covering fraction of 0.4 and pure 
thermal velocity for the lines. The determination of these pa¬ 
rameters must await better S/N observations. 

We started by re-analyzing the data obtained during the 
brightest flux state, that obtained from the ASCA observations in 
1994. There is a small flux change during the observations and 
our analysis refers only to the first 45 ks representing the high¬ 
est flux. The fit result is very similar to the one obtained by G98 
(model B(I) table 5 in their paper). It requires a single full line- 
of-sight covering HIG cloud with a column density of Nh = 
lO^'-'^ cm-2, Uox=l0-^-^ and the r(0.5 -50keV) = 1.9 SED. 
We assumed a neutral galactic column of 3x10^*^ cm^^(see 
G98) and constant gas density of «// = 10^ cm^^. As argued 
below, this is consistent with our lower limit estimate of the 
HIG density. None of the following conclusions depends much 
on the precise value of «//. Tests show that the uncertainty on 
the 0.5-50 keV spectral index is well below 0.2. 

The initial fit of the Chandra observations was based 
on the assumption of no change in SED and a large luminosity 
drop of the source as inferred from the hard X-ray observations. 
The basic premise is that all spectral variations are due entirely 
to flux variations that take place on a time scale much shorter 
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than the seven year period between observations. The fit ob¬ 
tained in this way compares poorly with the observations due 
to a large excess of flux around 5-7 keV not represented in the 
model. This can be explained by the fact that we also observed 
a strong, unresolved (o < 6600 km/sec) Ee 6.4 keV line with 
EW= 290 ±50 eV. This line intensity is clearly too strong to 
be emitted by the absorbing gas at its 2000 level of ionization 
and must therefore represent emission from another, unrelated 
component. The most likely explanation is that the line is due to 
emission by low ionization gas (i.e. much less ionized than the 
HIG), such as the one commonly observed by “reflection” in 
type-II (narrow emission lines) AGN. The typical Ee 6.4 keV 
EW in type-II sources is about I keV. The theoretical value 
(Matt, Brandt & Eabian 1996; Netzer, Turner & George 1998) 
is also 1 keV and depends on the gas composition, the contin¬ 
uum slope and the viewing angle of the reflector. Eor the cho¬ 
sen SED, we expect an EW of 700-1000 eV relative to a pure 
scattering continuum assuming solar metalicity and intermedi¬ 
ate Compton depth reflector. Given the observed Ee 6.4 keV 
line in NGC 3516, we suggest that its origin is the low ioniza¬ 
tion reflecting gas and the observed-vs-predicted EW suggests 
that reflection from this component contributed between 30 and 
50 percent to the continuum at 6.4 keV during the Chandra ob¬ 
servations. Thus, the factor 5 decrease in the continuum flux 
observed at 6.4 keV implies a decrease in the true luminosity 
of the central source by a factor 8 to 10 (if 30% or 50% of the 
observed continuum is due to reflection, respectively). 

To further test the idea of a low ionization component 
to the Ee 6.4 keV line, we calculated the emission and reflec¬ 
tion spectrum of gas with a large column density (Nh = 10^^ 
cm^^) and U„x=l0^^'^■ Since the location of this component 
is probably several light-years from the center, and the variabil¬ 
ity time scale of the central source is several weeks or less, we 
had to make an additional assumption about the time-averaged 
flux level of NGC 3516. The present model assumes this to 
be similar to the high 1994 flux level. Given this level, the lu¬ 
minosity of the reflection component requires a covering frac¬ 
tion of Ci/A% ~ 0.4 — 0.7; a reasonable fraction given the likely 
geometry. The emission and scattering from this gas must be 
included when fitting the spectrum at all other epochs. In par¬ 
ticular, its contribution is predicted to be non-negligible even at 
maximum luminosity. We have therefore included it in our final 
model of the 1994 spectrum and, while the overall fit results are 
unchanged in terms of SED and column density, we note that 
between 6 and 10 percent of the flux at around 6.4 keV is due 
to this component. The fit results are given in Table 2. 


Table 2: fit parameters 


Observation 


d.o.f. 

Uox 

r(0.5-50 keV) 

ASCA 94 (high) 

1.07 

432 

-2.1 

1.9 

ASCA 95 (high) 

1.06 

429 

-2.3 

1.9 

ASCA 98 (high state) 

0.96 

373 

-2.5 

1.7 

ASCA 98 (low state) 

0.68 

373 

-2.6 

1.7 

Chandra 0th order 

0.91 

263 

-3.2 

1.7 

Chandra ACIS/LETGS 

0.48 

1932 

-3.2 

1.7 


Given this additional reflection component, with scat¬ 
tered continuum and unresolved Ee 6.4 keV line, we have at¬ 
tempted another spectral fitting of the Chandra observations 
with the above SED. The fit is statistically acceptable (see Ta¬ 
ble 2) suggesting that the reflection component is, indeed, very 
important for explaining the spectrum at low flux levels. Eor 


reasons related to the 1998 ASCA observations (see below) we 
have also attempted a fit with the r(0.5 — 50 keV) = 1.7 SED. 
This continuum is also consistent with the 2000 Chandra ob¬ 
servations. Eig. 2 shows the r(0.5 — 50 keV) = 1.7 SED model 
on top of the new observations. Inspection of the diagram, and 
standard statistical analysis (Table 2) clearly demonstrate the 
applicability of this approach (note however the discrepancy 
over the 1-2 keV range are due to the known LETGS calibra¬ 
tion feature). In particular, we conclude that the very large flux 
drop, by a factor of ^ 50 below 1 keV between 1994 and 2000, 
is entirely consistent with the constant SED continuum and the 
dramatic increase in the absorption opacity due to the decreas¬ 
ing level of ionization. 

The noticeable flux variations correspond, according to 
the model, to a very large change in fractional ionization of the 
gas. While the strongest absorption features in 1994 are due to 
Ovil and Ovill, most of the opacity in the 0.6-0.9 keV energy 
range, in 2000, is due to Ov, Ovi and Ovil. According to the 
model, the continuum role-over below 0.6 keV is mostly due to 
Z < 8 elements, in particular carbon and helium. In addition, 
the model successfully predicts the observed N VI 0.419 keV 
and O VII 0.561 keV emission line intensities, given the obser¬ 
vational uncertainties and assuming the above mentioned cov¬ 
ering fraction by HIG clouds with properties similar to those of 
the line-of-sight material. This by itself is not a conclusive ev¬ 
idence for the absorber’s properties since the low ionization re¬ 
flecting component contribution to the O VII 0.561 keV line can 
approach the observed line intensity given the uncertain level 
of ionization of the reflecting gas. 

3.2. Spectral analysis of the 1995 and 1998 observations 

The seven-year history of ASCA, BeppoSAX and Chan¬ 
dra observations of NGC 3516 provide a unique opportunity 
to test the suggestion that all observed spectral variations are 
caused by a single variable, the unattenuated X-ray flux. To 
this end, we have collected most of the previous ASCA obser¬ 
vations of the galaxy and fitted, each one, in the same way as 
described for the 2000 data. The only data set not to be ana¬ 
lyzed is from December 1999. These data suffer from severe 
calibration problems at low energies where accurate spectra are 
required to test our model. 

The results of fitting the ASCA data are as follows; 

The 1995 data set: The mean 0.7-10 keV 1995 luminosity is 
about a factor 1.5 lower than the 1994 luminosity and 
the analysis refers to the first, higher flux part of the ob¬ 
servations to avoid the continuum variations. We have 
kept the r(0.5-50 keV)= 1.9 SED and scaled Uox by this 
factor to obtain a satisfactory fit (Table 2) to the data. 
Thus the observed 1995 spectrum is entirely consistent 
with our previous assumption. 

The 1998 data set: NGC 3516 was monitored continuously 
for about 360 ks during 1998 (see Nandra et al. 1999 
for details of this campaign). The average 6.4 keV con¬ 
tinuum flux was about 2.5 higher than the Chandra 2000 
flux, suggesting an intrinsic hard X-ray continuum some 
times more luminous than the 2000 continuum. The 
observations show a large (50%) peak-to-peak varia¬ 
tions on a time scale of several hours. 

Our attempt to fit the mean 1998 spectrum with the 
r(0.5 — 50 keV) = 1.9 SED and the same obscuring 
column, yielded unsatisfactory results. We found that. 
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depending on reflector luminosity, the E >2 keV slope 
is flatter than in 1994 and 1995, with 1.7 < E < 1.8. The 
result is unchanged when splitting the data into high and 
low states. We conclude that the source has undergone 
a real change in slope but by no more than AT = 0.2, 
between 1995 and 1998 (we re-emphasize that the 1994 
0.5-10 keV continuum is significantly steeper than F = 
1.7). However, assuming the r(0.5 —50 keV) = 1.7 
SED and normalizing Vox by the continuum flux ratio 
at 6.4 keV between 1998 and 2000, we can obtain satis¬ 
factory fits to the high and low 1998 spectra. 

The small change of slope at high energies in 1998 is the 
main motivation for considering two, slightly different 
SEDs for fitting the 2000 observations. The quality of 
the Chandra spectra is not good enough to distinguish 
between the two since both give satisfactory fit to the 
observations. However, we should note that the success 
of both models depend crucially on the assumed 0.1-0.5 
luminosity since the (observationally unknown) flux in 
this band determines the opacity through the ionization 
of the Z < 8 elements. The assumed 0.1-0.5 keV slope 
for both SEDs is a direct result of those constraints and 
the need to adequately fit the data above 0.5 keV. 

To summarize, our model fittings suggest r(0.5-50 
keV)=1.9 in 1994 and 1995 becoming flatter, by AT = 0.1—0.2 
in 1998. The 2000 observations are consistent with either one 
of those SEDs. Eig. 3 shows our best models and data-to-model 
ratios for all data sets (1994, mean-1995, low-1998 and 2000). 
Table 2 gives the fit parameters and the calculated all in 
agreement with our acceptability criteria. 

3.3. Color Analysis of the 1998 Data 

Given the satisfactory spectral fits, we can now look 
for the shortest variability time corresponding to a detectable 
change in the gas level of ionization. The 1998 data set is most 
suitable for this since the high energy continuum varied sub¬ 
stantially during the observations. Unfortunately, it is difficult 
to obtain meaningful constraints on the time scales involved us¬ 
ing complete spectral fitting since it requires integrating over 
at least 20-30 ks; a period which is long considering the ob¬ 
served variations in 1998. Therefore we adopted the following 
approach; we have fitted the spectra at two epochs represent¬ 
ing local minima and maxima separated by 60 ks (see details 
below) and we have also used a complementary method based 
on broad-band photometry. The broad-bands approach follows 
closely the “color analysis” method discussed in Netzer, Turner 
& George (1994) and shown to be useful for analyzing absorbed 
X-ray spectra. The method allows the analysis of much shorter 
data trains with much improved S/N and considerably better 
time resolution. We chose two bands, XMl which gives the in¬ 
tegrated flux over the 0.7-1.2 keV range and XM3 that covers 
the 4-10 keV range. The fluxes in these band were measured 
for the entire 1998 data set in bins corresponding to one orbit 
(~ 4000 sec on target integration) and the “color” (or “softness 
ratio”) XM1/XM3 was calculated for each of those . The results 
are presented as a function of time, in Eig. 4. 

Inspection of Eig. 4 demonstrates the close agreement 
of the XMl and XM3 flux variations. Cross correlation analy¬ 
sis shows the two bands to be highly correlated, with no mea¬ 
surable time lag, confirming the earlier result of Edelson et al. 
(2000). This tight correlation must be due to intrinsic contin¬ 
uum flux variations. The noticeable changes in color, over time 


scales of a few orbits, are also very clear. We suggest that they 
are caused by real changes in the absorber’s opacity because of 
the following reasons: Eirst, we have verified, by detailed spec¬ 
tral fitting, that the change of hard continuum spectral index be¬ 
tween minimum and maximum flux levels during 1998 is less 
than AT = 0.1. This corresponds to a too small change (factor 
of less than 1.1) to be responsible to the observed variations in 
XM1/XM3 (factor about 1.3). Second, our model calculations 
(see Eig. 4) are entirely consistent with the variability ampli¬ 
tude of XM1/XM3 confirming that the changing opacity is the 
main driver for the observed changes. Cross correlation analy¬ 
sis suggests a short delay between the XM3 and the XM1/XM3 
light curves but this is not statistically significant. 

We conclude that a significant change in level of ioniza¬ 
tion has occurred on a time scale of a few orbits. A conservative 
upper limit for the time of significant XM1/XM3 change is 60 
ks, which is roughly the time between the flux level near t=240 
ks and the minimum near t=300 ks. Direct spectral fittings of 
these two epochs, assuming Vox variations in direct proportion 
to the observed flux changes, are entirely consistent with this 
idea (see Table 2). They produce acceptable fits and confirm 
that pure continuum slope variations cannot explain the ob¬ 
served effect (note again the limitation of spectral fitting that 
requires integration over many orbits to achieve adequate S/N). 
In particular, the calculations show that the dominant stage of 
ionization is Ovil, with fractional ionization of about 0.8, and 
most of the low energy opacity variations are due to the increas¬ 
ing opacity of Ovi following the continuum flux variations. 

The determination of tovi, the time associated with the 
increase in the Ovi abundance until it becomes a significant 
opacity source (i.e. to exceed a fractional ionization of 0.1), 
requires full time-dependent solution of the ionization and re¬ 
combination equations. In particular the commonly used 1/e 
time scale, referred to as “recombination time”, is not necessar¬ 
ily a good measure under such conditions. However, recombi¬ 
nation onto Ovi provides a good estimate for tovi, given known 
ionization equilibrium solution and the fractionl ionization of 
Ovi and Ovil. Using these fractions, the known recombination 
rates at T ~ 3.5 x 10^ K (the temperature obtained from the 
model) and neglecting di-electronic and three-body recombina¬ 
tion that are not important at these densities and temperatures, 
we find 


ne 

where T$ is the electron temperature in units of 10^ K. This 
gives a lower limit of about 2.4 x 10® cm^^ for the recombin¬ 
ing gas density. Combined with the known X-ray luminosity 
during the 1998 observations, we derive an upper limit on the 
distance of the absorbing gas of about 6 x lO'V cm. Given 
the range in density and location, and noting the conservative 
limits on both, the preferred geometry of the absorbing material 
is that of a thin shell with AR/R « 1 

4. DISCUSSION 
4.1. Alternative models 

While the agreement between the ~7 years of observa¬ 
tions and the model presented here is satisfactory, we have nev¬ 
ertheless considered various alternative explanations. 

Eirst, we considered the possibility that all observed 
changes are due entirely to variations in continuum slope. Spec¬ 
tral fits show that this requires a slope of F = 1.13 in 2000 
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Fig. 3.— Left: Adopted models for the 1994 (high level), 1995 (high level), 1998 (low level at t=300 ks, see Fig. 4) and 2000 observations. All models assume 
Nh = 10^*'^ and the scaling of Uox is according to the intrinsic continuum luminosity (see text). Right: data-over-model ratios for the models shown on the left. 
Both ACIS/LETGS data (in magenta) and the zeroth-order spectrum are shown for the Chandra observations 


and the same Uox as derived for the 1994 observations (i.e. no 
change in opacity between the two epochs). This scenario cor¬ 
responds to a case where the absorbing gas is of very low den¬ 
sity, such that its level of ionization reflects the average lumi¬ 
nosity of the source (which is similar to the 1994 luminosity 
given the required Uox)- The major problem is that, at this ion¬ 
ization parameter, the emission lines produced by such gas are 
much stronger and more highly ionized compared with the lines 
observed by Chandra. Such lines would easily be detected by 
the ACIS/LETGS observations. 

A second possibility is that the absorbing gas is in¬ 
deed responding to the incident flux yet there are also large 
slope variations associated with the large changes in luminos¬ 
ity. Slope variations are claimed to take part in many AGN 
(Markowitz and Edelson 2001 and references therein) and are 
well known in galactic X-ray sources. It is therefore important 
to consider their effect on the 1994-2000 spectra of NGC 3516. 

As illustrated by the detailed analysis of the 1998 and 
2000 spectra, and as confirmed by various other tests, we found 
it practically impossible to identify a consistent model with 
large spectral index variations that fits the new observations. 
The main reason is the observational constraints on the hard X- 
ray flux. Since the hard X-ray continuum is directly observed, 
we can easily calculate the change in Uox associated with any 
assumed change in E. Since the SED and the level of ionization 
of the gas, during 1994, are known, we can easily calculate the 
effect such changes will have on the line-of-sight gas opacity. 
As explained, our model assumed no (or very small) change in 
slope and hence its predictions represent the smallest possible 
changes in level of ionization. Thus, the assumption of addi¬ 
tional continuum flattening is translated to a larger change in 
Uox and hence larger opacity at times of lower luminosity. The 
computations show that such an effect is so large that it will 
produce a low luminosity spectrum which is inconsistent with 
the Chandra observations. In fact, as demonstrated earlier, the 
largest change of 0.5-50 keV spectral index that is consistent 
with the single absorber model is about AF = 0.2. 

Another consistency check on the model is the com¬ 
parison between UV and X-ray observations. NGC 3516 was 
observed by HST, HUT, and FUSE on numerous occasions. 



time (10® sec) 


Fig. 4. — Binned (two-orbits) XMl and XM3 fluxes, and XM1/XM3, as 
functions of time for the 1998 ASCA observations. A smooth model was fitted 
to the XM3 observations (solid line) and the model predictions are shown on 
the XM1/XM3 diagram. 

some of these simultaneous with ASCA and Chandra observa¬ 
tions. While a detailed discussion of the UV-X-ray relationship 
is deferred to a later paper, we note that the overall ranges in 
UV and hard X-ray luminosities is very similar indicating small 
variations in the SED. 

Einally, combined luminosity, SED and column density 
variations may also be considered. The number of free parame¬ 
ters in such models is so large that additional constraints, from 
a large number of emission and absorption lines, are required 
to assess their reality. 

Several other papers discussed the X-ray spectrum of 
NGC 3516 and suggested alternative explanations for its spec¬ 
tral variations. Kriss et al. (1996b) suggested a range of 
ionization parameters and column densities. Kolman et al. 
(1993) suggest a connection with the BLR and Mathur et al. 
(1997) discuss the UV-X-ray absorption connection. A recent 
BeppoSAX-hased study, reported in Guainazzi et al. (2001), dis¬ 
cusses two epochs during 1996-1997, separated by 7 months 
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and showing large spectral variations similar to the ones dis¬ 
cussed here. The authors adopted a completely different ap¬ 
proach by considering both warm and cold absorbers with 
different column densities. Their interpretation of the two 
different spectra is very different from ours involving non¬ 
equilibrium situations, similar to the model proposed by Nicas- 
tro et al. (1999), and large changes of spectral index. Their fa¬ 
vorable model includes cold absorbing clouds of large columns 
{few X 10^^ cm^^) moving in and out of the line of sight. Ac¬ 
cording to the model, the observed variations are due to the 
onset of substantial cold absorption during times in which the 
properties of the primary continuum do not change. 

There are observational limitations and theoretical prob¬ 
lems in the Guainazzi et al. (2001) interpretation. The use 
of absorption edges instead of a self consistent photoionization 
model, and the assumption of a large change of slope, prevented 
them from testing the simpler idea that the apparent changes in 
spectral index are due to changes in the level of ionization of 
the absorber. Visual inspection of their data (Fig. 1 there) sug¬ 
gests, infact, that the two BeppoSAX observations can be fitted 
into the same sequence of ionization vs. luminosity presented 
in Fig. 3. 

4.2. Ultraviolet emission and absorption lines 

The model proposed here has important observational 
consequences to the UV properties of the source at minimum 
light. Given the model assumptions we find large column densi¬ 
ties of Ovi, Nv and CiV and hence strong emission and absorp¬ 
tion lines produced by these species during the Chandra obser¬ 
vations. The exact values depend on the assumed (unobserved) 
SED and will be discussed, together with new HST/STIS data, 
in a future publication. 

Mathur et al. (1997) discussed possible relationships be¬ 
tween X-ray and UV absorbers in NGC 3516. Their general 
picture is similar in some ways to the global properties of the 
source discussed here, but the actual values they have derived 
for the absorbing column and the ionization parameter are quite 
different from ours. Much of the difference is likely due to the 
limited resolution of the ROSAT-PSPC 1992 data they used for 
their analysis. They have noticed a large change of slope in ear¬ 
lier observations of the source but did not interpret it as due to 
the change in opacity. 

An interesting aspect of the Mathur et al. study is the 
suggested link between the EUV and X-ray properties of the 
absorber. In particular, they studied the possibility of an out¬ 
flowing absorbing gas that produce both UV and X-ray absorp¬ 
tion and made a clear prediction of the “disappearance” of such 
gas, over a period of ~ 15 years, due to the constant increase in 
its level of ionization. Our Chandra observations are in clear 
contradiction to this prediction suggesting that the increased 
opacity in 2000 may be from the same gas responsible for the 
strong absorption a decade earlier. 

4.3. The location, density and mass of the HIG 

Our results produce the best limits obtained so far on the 
density and location of the HIG in type-I AGN since they are 
based on the internal consistency of several observations fitted 
with a very simple model over a very large energy range. The 
only important variable of the model is the ionization param¬ 
eter which is normalized by the observed high energy contin¬ 
uum flux in 1994. As argued above, small intrinsic spectral 
index variations are probably seen but their effect is far smaller 


than the change due to absorber’s opacity. Our observations and 
model fittings show the absorber to lie well inside the central 
Ipc and perhaps much closer in, at the BLR location or perhaps 
even closer. 

There were several previous attempts to deduce the HIG 
location in AGN. Guainazzi et al. (1996) used a long-look 
ASCA observation of NGC 4051 to claim a significant change 
in the Ovil edge, over a time scale of 10^ seconds during 
which there was no change in the Ovill edge strength. How¬ 
ever, they did not show a realistic model to explain this change. 
Guainazzi et al. (1998) and Uttley et al. (1999) observed NGC 
4051 in a very low state that allowed them to place a lower 
limit on the distance of the cold reflector. Komossa & Meer- 
schweinchen (2000) investigated HIG systems in four radio- 
loud quasars and NFS Is, using ROSAT data, and suggested 
edge variations. While the claimed variability is probably real, 
the ROSAT-PSPC resolution does not allow accurate modeling 
of the absorption features. Moreover, so far there is no posi¬ 
tive detection of HIG with “typical” AGN absorption features in 
any NLSl. Reynolds and Fabian (1995) and Otani et al. (1996) 
made detailed studies of the absorption features in MCG-6-30- 
15 and claimed an Ovill edge change on a time scale of 10“* 
sec. They also claimed that the Ovil edge was observed but 
did not change during the same time which caused them to pro¬ 
pose that this feature is due to gas much further away from the 
central continuum. While the observed variations and the as¬ 
sociated time scales are possible, the interpretation is highly 
questionable in view of the new Chandra and XMM-Newton 
observations of this source (Branduardi et al. 2000; Lee et al. 
2001). In particular, it is not at all clear whether the observed 
0.7-0.9 keV absorption-like features are due to the Ovill and 
Ovil edges suggested by Otani et al.. Infact, they may be re¬ 
lated to a complex of broad emission lines and much narrower 
and weaker absorption lines (Branduardi et al. 2000). A more 
consistent physical model is certainly needed to understand this 
source. 

Our new findings are of great importance for AGN study 
since they are the first to show clear and direct changes of the 
absorber properties in response to changes in source luminos¬ 
ity. The time span is very long (7 years) and the variability 
amplitude very large (a factor of ^ 8). The lack of clear dis¬ 
tance indicators made it impossible, until now, to distinguish 
between HIG systems observed in absorption in type-I sources 
and the ones observed in emission in type-II sources. Given the 
deduced location in NGC 3516, and a plausible geometry of the 
obscuring medium (torus?) in type-II sources, it is evident that 
the system we have observed is very different, in location and 
other properties, from the extended X-ray emission line gas in 
type-II AGN. Unfortunately, there is no precise mass estimate 
for this gas. The combined density and location we derived 
give an upper limit to the mass of about 80 M© (which is or¬ 
ders of magnitudes larger than the mass of the ionized gas in 
the BLR) but much smaller values are entirely consistent with 
this estimate. 

4.4. Neutral reflectors and relativistic Fe 6.4 keV lines 

Important ingredients in our model are the line and con¬ 
tinuum contributions from the low ionization reflector. As ex¬ 
plained, this component is necessary to explain the observed 
Fe 6.4 keV at minimum flux and is consistent with the idea of 
a large column-density torus like the ones suggested for type- 
II sources. Interestingly, this component contributes a signifi- 
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cant amount to the observed flux even at maximum (1994 level) 
luminosity: about 10% of the 6.4 keV flux and EW(narrow 
Ee 6.4 keV)~ 30 eV under the assumption of 50% contribution 
at 6.4 keV during the Chandra observations. This has impor¬ 
tant consequences to the intensity and profile of the so-called 
“relativistic” Ee 6.4 keV line in NGC 3516 and other AGN. In 
particular, the claimed Ee Ka line in this source (e.g. Nandra et 
al. 1999) has a typical broad red shoulder and a sharp blue side. 
The cold gas contribution, if similar to that found here, peaks 
at 5-1 keV and will have a significant effect on the relativistic 
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